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Rectified transport of active ellipsoidal particles is numerically investigated in a two-dimensional 
asymmetric potential. The out-of-equilibrium condition for the active particle is an intrinsic prop¬ 
erty, which can break thermodynamical equilibrium and induce the directed transport. It is found 
that the perfect sphere particle can facilitate the rectification, while the needlelike particle destroys 
the directed transport. There exist optimized values of the parameters (the self-propelled veloc¬ 
ity, the torque acting on the body) at which the average velocity takes its maximal value. For 
the ellipsoidal particle with not large asymmetric parameter, the average velocity decreases with 
increasing the rotational diffusion rate, while for the needlelike particle (very large asymmetric pa¬ 
rameter), the average velocity is a peaked function of the rotational diffusion rate. By introducing 
a finite load, particles with different shapes (or different self-propelled velocities) will move to the 
opposite directions, which is able to separate particles of different shapes (or different self-propelled 
velocities). 
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I. INTRODUCTION 

Noise-induced transport far from equilibrium plays a 
crucial role in many processes from physical and bi¬ 
ological to social systems. The transport properties 
of systems consisting of active particles have generated 
much attention. There are numerous realizations of ac¬ 
tive particles^ in nature ranging from bacterial— and 
spermatozoa^ to artificial colloidal microswimmers. Self¬ 
propulsion is an essential feature of most living systems, 
which can maintain metabolism and perform movement. 
The kinetic of self-propulsion particles moving in poten¬ 
tials could exhibit peculiar behavior— The problem 
of rectifying motion in random environments is an im¬ 
portant issue, which has many theoretical and practi¬ 
cal implications^. At equilibrium the periodic potential 
alone is not able to produce a rectification effect, due to 
the detailed balance preventing time symmetry breaking. 
Indeed, one has to add some perturbation which breaks 
the time symmetry and brings the system out of equi¬ 
librium. For the active particles, the out-of-equilibrium 
condition is an intrinsic property of the system. So the 
active particle without any external forces can break the 
symmetry of the system and be rectified in periodic sys¬ 
tems. 

Recently, rectification of self-propelled particles in 
asymmetric external potentials has attracted much at¬ 
tention. Angelani and co-workers^i studied the run-and 
tumble particles in periodic potentials and found that the 
asymmetric potential produces a net drift speed. Even 
in the symmetric potential a spatially modulated self¬ 
propulsion and a phase shift against the potential can 
induce the directed transport^. Recently, transport of 
Janus particles in periodically compartmentalized chan¬ 
nel is investigated^ and the rectification can be orders 


of magnitude stronger than that for ordinary thermal 
potential ratchets. In all these studies^ - — on active 
ratchet, the active particle was treated as the point spher¬ 
ical particle. 

However, shape deformation of particles plays an im¬ 
portant role in nonequilibrium transport processes2Zr~. 
Han and co-workers2I experimentally studied the Brow¬ 
nian motion of isolated ellipsoid particles in two di¬ 
mensions and quantified the crossover from short-time 
anisotropic to long-time isotropic diffusion. In the pres¬ 
ence of an external potential, the external force can am¬ 
plify the non-Gaussian character of the spatial probabil¬ 
ity distributions^. Ohta and co-workers^ found that an 
isolated deformable particle exhibits a bifurcation such 
that a straight motion becomes unstable and a circu¬ 
lar motion appears. Due to the coupling of the ro¬ 
tational and translational motion, Brownian motion of 
asymmetrical particles is considerably more complicated 
compared to the spherical case, and thus shows peculiar 
behavior. Therefore, how active asymmetrical particles 
are rectified from a ratchet potential may receive much 
attention. In this paper, we will extend the study of ac¬ 
tive ratchet from the spherical particle to the asymmetric 
particle. We focus on finding how the asymmetry of the 
particle affects the rectified transport and how asymmet¬ 
ric particles can be separated. 


II. MODEL AND METHODS 

We consider an ellipsoidal (anisotropic) particle mov¬ 
ing in a two-dimensional ratchet potential (shown in Fig. 
1). The particle is self-propelled along its long axis. In 
the lab x-y frame, the particle at a given time t can be de¬ 
scribed by the position vector R(t) of its center of mass, 
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where vq is the self-propelled velocity in the body frame, 
which is taken along the long axis of the ellipsoid particle. 
Ti and T 2 are the mobilities along its long and short axis, 
respectively. T 3 is the rotational mobility and r is the 
torque acting on the body due to it orientation relative 
to the direction of the potential. F x and F v are the forces 
along x and y direction of the lab frame. The noise £ ?; (t) 
has mean zero and satisfies 

(£i{t)£j(t )) = - t ),i,j = 1,2,3, (4) 

1 i 

where T is the temperature and kg is the Boltzmann 
constant. 

We now obtain these equations in the lab frame based 
on the method described in Reb^S. By means of a 
straight forward rotation of coordinates, the displace¬ 
ment in the two frames are related by the following equa¬ 
tions, 


Sx = cos 98x — sin 9Sy, 


( 5 ) 


Sy = sin 98x + cos 98y. 


( 6 ) 


After somewhat manipulations, Eqs. mm in the 
body frame can be replaced by the following equations 
in the lab fram o 27 ' 28 


dx 

— = v 0 cos 9(t)+F x 
dt 


T+AT cos 20(f) 


+ATF y sin20(f)+£i(f), 

( 7 ) 


FIG. 1: Schematic diagram illustrating the system, (a) The 
active particle with the self-propelled velocity moving in the 
ratchet potential and the equipotentials of the potential de¬ 
scribed by Eci. (li3l) . (b)Representation of an ellipsoid in the 
x-y lab frame and the x-y body frame. The angle between 
two frames is 9. The displacement R can be decomposed as 
(Sx, Sy) or (Sx, Sy). 


dy_ 

dt 


vo sin 9(t)+F y 


T—AT cos29(t) 


+A TF X s'm29(t)+^2(t), 


( 8 ) 


d d(t) 
dt 


r3T + ^3 (t), 


( 9 ) 


which also corresponds to the coordinates in the body 
frame (Sx, Sy). 8(t) is the angle between the x axis of 
the lab frame and the x axis of the body frame. Ro¬ 
tational and translational motion in the body frame are 
always decoupled, so the dynamics of the active ellipsoid 
particle is described by the Langevin equations in this 
fram o 27 ' 28 , 


where the quantities f = ^(^-1-^) and AT = ^(Ti—T 2 ) 
are the average and difference mobilities of the body, re¬ 
spectively. The parameter AT determines the asymmetry 
of the body, the particle is a perfect sphere for AT = 0 
and a very needlelike ellipsoid for AT —» T. The noise 
£i(t) has mean zero and the following relation s 27 ' 28 

(&(fi)&(t))=2D e 6{t-i), (10) 


1 dx 

T\~dt 


= F x cos 9(t) + F y sin 9(t) + (t) + 
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(1) 


=2k B TT l3 S(t-t), 


J _dy_ 
r 2 dt 


F y cos 9(t) - F x sin 9(t) + | 2 (t), 


and 

( 2 ) 


Tij — T Sij + AT 


cos 29 sin 29 \ 
sin 29 — cos 29 J ’ 


( 11 ) 


( 12 ) 


1 d 9(t) 

r \~dT 


T + £z{t), 


where Dg = ksTT^ is the rotational diffusion rate, which 
describes the nonequilibrium angular fluctuation. The 
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statistical averages have superscripts to indicate over 
which noise is the average taken and subscripts to de¬ 
note quantities which are kept fixed. 

For the asymmetric potential, we consider the follow¬ 
ing potential^ (shown in Fig. 1(a)), 


U(x, y) = -^-y 2 [ cos(x + A In cosh y) + 1.1] + fx, (13) 

where Uq is the height of the potential and / is the load 
along the x direction. A is the asymmetric parameter of 
the potential and the potential is symmetric at A = 0.0. 
The equipotentials are now symmetry broken and look 
like a herringbone pattern for A ^ 0. 

In this paper, we focus on the direction transport 
of active asymmetrical particles. The behavior of the 
quantities of interest can be corroborated by Brown¬ 
ian dynamic simulations performed by integration of 
the Langevin equations mm using the second-order 
stochastic Runge-Kutta algorithm. Because the parti¬ 
cle along the y-direction is confined, we only calculate 
the x-direction average velocity based on Eqs. ms®, 


V9 0 


= lim ■ 

t—> OO 


&(*)> 


flifa 


(14) 


where 9q is initial angle of the trajectory. The full average 
velocity after a second average over all 9 0 is 

1 f 2 ^ 

V = —J d9 0 vg 0 . (15) 

For the convenience of discussion, we define the scaled 
average velocity v s = v/vo through the paper. In our 
simulations, the integration step time At was chosen to 
be smaller than 10~ 4 and the total integration time was 
more than 3 x 10 5 and the transient effects were esti¬ 
mated and subtracted. The stochastic averages reported 
above were obtained as ensemble averages over 3 x 10 4 
trajectories with random initial conditions. 
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FIG. 2: Scaled average velocity v s as a function of the asym¬ 
metrical parameter A of the potential for different values of 
Ar at vo = 2.0 and Dg = 0.1. 
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FIG. 3: Scaled average velocity v s as a function of the asym¬ 
metrical parameter Ar of the particle for different values of 
vo at A = 1.0, / = 0.0, and Dg = 0.1. 


III. RESULTS AND DISCUSSION 

Based on the numerical simulations, we mainly calcu¬ 
late the average velocity for the two cases: zero load and 
finite load. For the zero load case, we focus on the recti¬ 
fication effects and how the parameters can affect recti¬ 
fication. For the finite load case, we present two particle 
separation methods: shape separation and self-propelled 
velocity separation. In the simulations, unless otherwise 
noted, we set ksT = 1.0, r = 0.0, f = 1.0, and Uq = 1.0 
throughout the paper. 

A. Zero load and rectification 

We first consider the zero load case (/ = 0.0). The 
scaled average velocity v s as a function of asymmetrical 


parameter A of the potential is reported in Fig. 2. It 
is found that v s is positive for A > 0, zero at A = 0, 
and negative for A < 0. A qualitative explanation of 
this behavior can be given by the following argument. 
For A = 0 (symmetric case ) the probabilities of cross¬ 
ing right and left barriers are the same and then there is 
a null net particles flow. For A > 0, the left side from 
the minima of the potential is steeper, it is easier for 
particles to move toward the gentler slope side than to¬ 
ward the steeper side, so the average velocity is positive. 
Therefore, the asymmetry of the potential will determine 
the direction of the transport and no directed transport 
occurs in a symmetric potential. 

The dependence of the scaled average velocity v s on 
the asymmetrical parameter AT of the particle is pre¬ 
sented in Fig. 3 at A = 1.0. We find that v s decreases 
monotonically with the increase of the asymmetrical pa¬ 
rameter AT. In order to give the explanation of the 
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FIG. 4: Scaled average velocity v 3 as a function of the self- 
propelled velocity vo for different values of Ar at A = 1.0,/ = 
0.0, and Dg = 0.1. 


phenomenon, we present the translational diffusion co- 
efficient in the x direction^ D xx = D 0 + —— -" . 
where Do = ksTT. As we know, the increase of D xx en¬ 
hances the ratchet effect when D xx < Uo and reduces the 
ratchet effect when D xx > Uo- When D xx ~ Uq, the opti¬ 
mized ratchet effect occurs. In our system, Dq = 1.0 and 
Uo = 1-0, the ratchet effect is optimized when Ar = 0 
(D xx = U 0 ). As Ar increases from zero, D xx > Uq, 
the ratchet effect is gradually destroyed and v s decreases 
monotonically. Therefore, the perfect sphere particle can 
facilitate the rectification, while the needlelike ellipsoid 
particle destroys the directed transport. 

From Fig. 3, we can also find that the curve is con¬ 
vex for large value no and concave for small value vo- 
This phenomenon can be easily explained by introduc¬ 
ing the two factors: (A) the increase of no (from 0 to 2) 
enhances the transport (shown in Fig. 4) and (B) the 
increase of AT reduces the transport. For small value 
n 0 (no = 1.0), on increasing AT, factor B firstly domi¬ 
nates the transport, v s reduces quickly, and finally the 
ratchet effect gradually disappears (very small value v s ), 
v s reduces slowly, so the curve is concave. For no = 2.0, 
on increasing AT, factor A firstly determines the trans¬ 
port, v s reduces slowly, and finally factor B also becomes 
important, v s reduces quickly, so the curve is convex. 

Figure 4 shows the scaled average velocity n s as a func¬ 
tion of the self-propelled velocity no- The term no cos 6(t) 
in Eq. 0 can be seen as the external driving force. 
When no —> 0, the external driving force can be negligi¬ 
ble, so the average velocity will tend to zero. For very 
large values of no, the effect of the asymmetry of the po¬ 
tential reduces, thus v s becomes small. Therefore, there 
exists an optimal value of no at which v s takes its max¬ 
imal value. So the optimal self-propelled velocity can 
facilitate the rectification of particles. 

Results for v s as a function of the rotational diffu- 



FIG. 5: Scaled average velocity v s as a function of the 
rotational diffusion rate Dg for different values of Ar. 
(a)ellipsoidal particles (Ar = 0,0.3, 0.6). (b)needlelike parti¬ 
cles (AF = 0.9, 0.95, 0.99). The other parameters are A = 1.0, 
/ = 0.0, and no = 2.0. 


sion rate Dg are presented in Fig. 5 for different val¬ 
ues of Ar. For the case of ellipsoidal particles ( AT = 
0,0.3, 0.6) shown in Fig. 5(a), the scaled average ve¬ 
locity v s decreases with the increase of Dg. The term 
no cos 9(t)(cx no cos (Dgt)) in Eq. 0 can be seen as the 
external driving force along ^-direction. In the adiabatic 
limit Dg —> 0, the external force can be expressed by 
two opposite static forces no and —no, yielding the mean 
velocity V = |[n(no) + n(—no)], which is similar to the 
adiabatic case in the forced thermal ratchet^. As Dg in¬ 
creases, the scaled average velocity v s decreases. When 
Dg —>• oo, the self-propelled angle changes very fast, par¬ 
ticles are trapped in the valley of the potential, so n s 
tends to zero, which is similar to the high frequency driv¬ 
ing case in the forced thermal ratchet^. 

However, for the case of needlelike particles (AT = 
0.9,0.95,0.99) shown in Fig 5 (b), the anisotropic diffu¬ 
sion dominates the transport and there exists an optimal 
value of Ar at which v s takes its maximal value. This 
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FIG. 6: Scaled average velocity v a as a function of the torque 
r acting on the body for different values of Vo at A = 1.0,/ = 
0.0, Dg = 0.1, and AF = 0.0. 

is due to the mutual interplay between the anisotropic 
diffusion and the rotational diffusion rate. There are two 
time periods in the system: the anisotropic diffusion time 
T£> = jy— for crossing one period of the potential along 
x-direction and the period rg = A- for direction ran¬ 
domly varying in time. For the case of needlelike parti¬ 
cles, the anisotropic diffusion time becomes very impor¬ 
tant. When these two time periods cooperate with each 
other, the optimized rectification will occur. 

Figure 6 describes the dependence of the scaled aver¬ 
age velocity v s on the torque r for different values of vo- 
From Eq. ED, we can see that the self-propelled angle 
d(t) is determined by the torque and the random noise. 
When r —>■ 0, the rotational diffusion rate determines 
the angle 9(t). On increasing r, both the torque and the 
random noise play the important roles and work together 
in the transport, which induces the maximal average ve¬ 
locity. However, when r —> oo, the self-propelled angle 
changes very fast, particles are trapped in the valley of 
the potential and v s goes to zero. Therefore, large torque 
would suppress strongly the ratchet effect. 

B. Finite load and particle separation 

Since transport behaviors in the present system 
strongly depend on the asymmetric parameter AT of the 
particle and the self-propelled velocity vq , it is possible 
to realize particle separation. We will present two par¬ 
ticle separation mechanisms which induce the motion of 
particles of different AT or vo in opposite directions by 
introducing an external load / on the x-direction. 

Figure 7 (a) shows the scaled average velocity as a 
function of the load / for different values of AT. For 
small values of /, the positive ratchet effect dominates 
the transport and the average velocity is positive. On 
increasing the load /, the load dominates the transport, 



FIG. 7: (a) Scaled average velocity v a as a function of the 
load / for different values of Ar. (b)Scaled average velocity 
v a as a function of the parameter Ar at / = 0.9. The other 
parameters are A = 1.0 and vo = 2.0, and Dg = 0.1. 


the average velocity crosses zero and subsequently re¬ 
verses its direction. For very large load, particles are 
blocked and v s tends to zero exponentially. From Fig. 
1(a), we can find that if the load / is negative (positive 
force along x-direction), the particle will move forward 
without problem. If the load / is positive (negative force 
along x-direction), the particle will move backward and 
may get into a spine of the herringbone. Thus, the parti¬ 
cle has to go against the force in order to climb back up 
and get again on the backbone. Moreover, one can find 
that the stronger the force, the deeper the spine. There¬ 
fore, the particle is blocked in the spines for very large 
values of /. Note that this blocked phenomenon had 
been explained detailedly in Ref. 33 . There exists a valley 
in the curve v s — / at which the average velocity takes 
its negative maximal value. The position of the valley 
varies with AT. For a given load (/ = 0.9) shown in Fig. 
7 (b), the asymmetric parameter Ar of the particle can 
determine the direction of the transport. Particles larger 
than the threshold asymmetric parameter AT C move to 
the left, whereas particles smaller than that move to the 
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FIG. 8: (a) Scaled average velocity v s as a function of the 
load / for different values of vo- (b) Average velocity v as a 
function of vo at / = 0.75. The other parameters are A = 1.0 
and AF = 0.5, and Dg = 0.1. 


right. Therefore, one can separate particles of different 
values of AT and make them move in opposite directions. 

Figure 8 (a) shows the scaled average velocity as a func¬ 
tion of the load / for different values of the self-propelled 
velocity vq. There exist two driving factors in the sys¬ 
tem: (1) the self-propelled velocity vq which induces the 
positive current; (2) the load / which causes the neg¬ 
ative current. For small values of vq (e. g. vq = 0.5), 
the load / dominates the transport and the self-propelled 
velocity vq can be neglected. When / —»• 0, the average 
velocity goes to zero. When increasing /, the amplitude 
of the negative velocity increases. However, for very large 
values of /, the particle is blocked and the average ve¬ 
locity tends to zero. Therefore, there exists a valley in 
the curve v s — f at which the average velocity takes its 
negative maximal value. However, on increasing vq, the 
self-propelled velocity driving factor becomes more im¬ 
portant and the positive current gradually dominates the 
average velocity. Therefore, the valley in the curve grad¬ 


ually disappears. Therefore, the self-propelled velocity 
vo strongly affects the transport and even determines its 
direction. For a given value of / (e. g. 0.75) shown 
in Fig. 8 (b), particles with larger than the threshold 
velocity v c move to the right, whereas particles smaller 
than that move to the left. Therefore, one can separate 
particles of different values of vq and make them move in 
opposite directions. 

Note that there are some other methods which can 
separate particles for our system. Particles can be sep¬ 
arated if their drift speeds are different, no need to 
move in the opposite directions. Particles with differ¬ 
ent shapes can be separated by many other means rather 
than the ratchet method. However, our separation meth¬ 
ods (which make particles move in opposite directions) 
are more effective and feasible. 


IV. CONCLUDING REMARKS 


In this work we have studied the transport of active 
ellipsoidal particles in a two-dimensional asymmetric po¬ 
tential by numerical simulations. It is found that the 
self-propelled velocity can break thermodynamical equi¬ 
librium and induce directed transport. The direction of 
the transport is determined by the asymmetry of the po¬ 
tential. The shape of particles can strongly affect the 
rectified transport, the perfect sphere particle can facil¬ 
itate the rectification, while the needlelike particle de¬ 
stroys the directed transport. There exist optimized val¬ 
ues of the parameters (the self-propelled velocity, the 
torque acting on the body) at which the average velocity 
takes its maximal value. For the ellipsoidal particle with 
not large asymmetric parameter, the average velocity de¬ 
creases with increasing the rotational diffusion rate, while 
for the needlelike particle (very large AT ), the average 
velocity is a peaked function of the rotational diffusion 
rate. In addition, by introducing a finite load on parti¬ 
cles, we have presented two particle separation ways (1) 
shape separation: particles larger than the critical asym¬ 
metric parameter AT C move to the left, whereas particles 
smaller than that move to the right; (2)self-propelled ve¬ 
locity separation: particles with larger than the threshold 
velocity v c move to the right, whereas particles smaller 
than that move to the left. Therefore, one can separate 
particles of different shapes (or different self-propelled 
velocities) and make them move in opposite directions. 
The results we have presented have a wide application in 
many systems, such as spontaneously moving oil droplets, 
motor proteins, bacterial swimmers, and motile cells. 
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